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Witnessing the active assembly phase of massive galaxies 
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ABSTRACT 

We present an analysis of ^60 000 massive (stellar mass M* > IO^Mq) galaxies out 
to z = 1 drawn from 55.2 deg 2 of the United Kingdom Infrared Telescope (UKIRT) 
Infrared Deep Sky Survey (UKIDSS) and the Sloan Digital Sky Survey (SDSS) II 
Supernova Survey. This is by far the largest survey of massive galaxies with robust 
mass estimates, based on infrared (if-band) photometry, reaching to the Universe at 
about half its present age. We find that the most massive (M+ > 10 Mq) galaxies 
have experienced rapid growth in number since z = 1, while the number densities 
of the less massive systems show rather mild evolution. Such a hierarchical trend 
of evolution is consistent with the predictions of the current semi-analytic galaxy 
formation model based on ACDM theory. While the majority of massive galaxies are 
red-sequence populations, we find that a considerable fraction of galaxies are blue star- 
forming galaxies. The blue fraction is smaller in more massive systems and decreases 
toward the local Universe, leaving the red, most massive galaxies at low redshifts, 
which would support the idea of active 'bottom-up' formation of these populations 
during < z < 1. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: 
formation - galaxies: luminosity function, mass function - galaxies: stellar content - 
cosmology: observations. 



1 INTRODUCTION 

Understanding the origin and evolution of galaxies, in par- 
ticular the most massive, is one of the major challenges in 
modern astrophysics. Many massive galaxies today are gi- 
ant early- type systems; hence the formation of spheroids 
should proceed to a certain extent in locked step with 
the mass assembly. The compelling theory of hierarchi- 
cal galaxy formation predicts that galaxies are assembled 
through successive mergers of smaller systems in overden- 
sities, or haloes, of h ypothetical cold dark matter (CDM) 
(|White fc Reeslll97Sl ). Massive galaxies therefore emerge in 
the last phase of the formation history. Alternatively, mas- 
sive galaxies could form through the rapid collapse of gas 
followed by a single prominent starburst at h i gh re dshifts 
l|Eggen. Lvnden-Bell. fc Sandagel 1 19621 ; iLarsonl 1 19751 ) . This 
'monolithic' scenario is supported by, for example, the tight 
colour-luminosity relation of early-type galaxies found in 
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galax y clusters (e.g.. iBower. Lucev. fc Elliil 19921 ; lEli is et ajj 
1 19971 ). 

While different evolution in different models makes dis- 
tant massive galaxies a unique test-bed for galaxy formation 
scenarios, observations have not yet provided evidence for 
the evolutionary path of those galaxies. The major obstacle 
in observations originates from the scarcity of galaxies at the 
high end of the galaxy mass function; it means that not only 
it is hard to find the population but also cosmic variance, 
the field-to-field variation of observed volume density arising 
from large-scale structure, is significant. In the last decade, 
many large programmes of optical-band imaging have been 
carried out, providing excellent data sets with which to 
investigate distant red old g alaxies in wide fields of sky 



exceeding a whole deg 2 (e.g., Bell et alJ 20041; Borch et all 



20061; ICimatti. Daddi. fc Rcnzini 2006; Will mer et ai1l2006l ; 



Brown et al.ll2007l ; lFaber et al-lbOOTl hThev consistently sug- 



gest that the total stellar mass locked in red galaxies with 
luminosities around L ~ L* , where L* is the characteris- 
tic luminosity of the luminosity function, has at least dou- 
bled since z ~ 1. Some of them also claim little growth in 
the number of very luminous galaxies well above L ~ V . 
However, while luminous red galaxies roughly correspond to 
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massive galaxies, it is not clear how well the evolution in the 
number of galaxies at the steep high end of the mass function 
is understood from these results, since the much more nu- 
merous, less massive galaxies with mass-to-luminosity ratios 
slightly less than average could easily dominate the observed 
numbers of luminous galaxies. The above authors also re- 
veal that a field of view of the order of a whole deg 2 is 
still not sufficient to conquer the uncertainty arising from 
cosmic variance for the high-end populations of the galaxy 
mass function. 

The advent of the United Kingdom Infrared Tele- 
scope (UKIRT) In fraredDeep Sky Survey (UKIDSS; 
lLawrence et al.ll2007l ) provides a unique opportunity to pro- 
duce an ideal sample to trace various aspects of massive 
galaxies in the distant Universe. Here we report the results 
of a A'-band survey with optical (u, g, r, i, z band) and near- 
infrared (Y, J, H band) photometry and optical spectra, 
focusing on massive (M* > 10 11 Mq) galaxies out to z = 1 
in an unprecedented large area covering 55.2 deg 2 . The K- 
band photo metry provides robust estimates of galaxy stellar 
masses (e.g-. lMatsuoka et al.ll2008h while the very large field 
of view significantly suppresses cosmic variance, which allow 
us to conduct a unique analysis of the mean properties of 
distant massive galaxies. 

This paper is organized as follows. In Section 2 we de- 
scribe the data sources and reduction process to extract the 
galaxy sample from the available data. Photometric redshifts 
and stellar masses are measured for each galaxy in Section 3. 
In Section 4, the number- density evolution of massive galax- 
ies and the associated uncertainties are explored. We then 
discuss the star-forming properties of galaxies and the com- 
patibility of the present results with previous measurements 
in Section 5. A summary follows in Section 6. Throughout 
this paper, we adopt the concordance cosmology of Ho = 70 
km s^ 1 Mpc -1 , Qm = 0.3, and Q,a = 0.7. Magnitudes are 
expressed in the Vega magnitude system for the UKIDSS 
near-infrared bands and in the AB magnitude system for 
the optical bands. 



2 DATA SOURCES AND REDUCTION 
2.1 Near-Infrared Photometry 

We extract from the Data Release 3 (DR3; Warren et al., in 
preparation) of the UKIDSS/Large Area Survey (LAS) the 
7f-band sources with right ascensions fro m l h 15 m to 3 h 6 m 
on the Sloan Digital Sky Survey fSDSS: lYork et al.ll2000h 
southern equatorial stripe (see Section 2.2). The range in 
right ascension is chosen so that the 7f-band observations 
are fairly complete within the sample area. More than half 
of the JC-band sources have also been observed in the Y, 
J and/or H bands. We exclude the sources assigned with 
serious quality flags corresponding to ppErrBits attributes 
larger than 31, or near (within 20 arcsec of) the detector 
edges of any exposure. We also exclude those sources near 
bright sources. This is achieved by searching for bright point 
and exte nded sources in the Two-Micron All-Sky Survey 
(2MASS; ISkrutskie et~ai1l2006l ) catalogues and rejecting all 
the LAS sources within sufficiently large distances of the 
bright 2MASS sources. The total effective area of the obser- 
vations defining our sample is 55.2 deg 2 . 



We retrieved all the images in our sample area from the 
DR3 data base, and used the SOURCE EXTRACTOR, version 2.5 
(Bcrtin & Arnouts! fl996) for magnitude measurements. The 
total magnitudes (m to t) of the sources are measured with the 
SOURCE EXTRACTOR total magnitude algorithm MAG_AUT0. We 
measure the aperture magnitudes (m ap ) with circular aper- 
tures of several diameters, 4.5, 3.3, 2.8 and 2.6 arcsec, which 
correspond to 20 kpc at z = 0.3, 0.5, 0.7 and 0.9, respec- 
tively. Since the seeing condition is generally superior in the 
UKIDSS LAS (^0.8 arcsec) to that in the optical obser- 
vations by the SDSS (~1.5 arcsec), the near-infrared im- 
ages were smoothed with the Gaussian kernel in such a way 
that the resultant full widths at half-maximum (FWHMs) 
of stellar profiles are similar to those in the optical images 
of the same field. The seeing measurements and smoothing 
were performed in each of the small rectangular subareas 
of approximately 9 x 13 arcsec 2 . Then we ran the SOURCE 
EXTRACTOR on the smoothed images to measure the aperture 
magnitudes of the sources. 

The detection completeness of the isT-band sources can 
be estimated by comparing the numbers of the LAS detec- 
tions with those of th e much deeper UKIDS S Deep Extra- 
galactic Survey (DXS; lLawrence et al.ll2007h in an overlap- 
ping field. This 0.6-deg 2 field is a part of the DXS VIMOS 4 
field, which is centred at RA 22 h 17 m , Dec. +00° 24' on the 
celestial equator. While the field is outside the RA range of 
our LAS fields, we confirmed that the evaluation sample of 
this field has a similar distribution of magnitudes and their 
errors to those for our actual sample. Below we show that 
the derived detection-completeness function reproduces the 
galaxy number counts from our sample, in excellent agree- 
ment with previous measurements, while it gives the lower 
limit of the detection completeness for our massive galaxies. 
We define our sample as being brighter than the limiting 
total magnitude K to t = 17.9 mag, where the detection com- 
pleteness is higher than 0.5. 

2.2 Optical Photometry 

We use the optical u, g, r, i and z-band images on the SDSS 
southern equatorial stripe. The stripe has be en observed re- 
peate dly in the SDSS-II Supernova Survey l|Frieman et al.l 
l200Sf ) during 2005- 2007, as well as in the original SDSS. 
We retrieved all the available images taken on the stripe 
from the SDSS Data Rel ease 6 Supplemental and Super- 
nova Survey data bases (I Adelman-McCarthv et ail I2008T ) ■ 
and stacked them in each of the five bands. The images 
observed in runs 2738 and 3325 are set apart from others, 
since they are taken with the standard survey conditions of 
the original SDSS and can work as the reference frames for 
stellar photometry. We measure, for each retrieved image, 
the mean and root-mean-square (rms) of the sky counts and 
the sky transparency at the observation by comparing the 
stellar photometry of the relevant frame with that of the ref- 
erence frames. After discarding the worst 5 per cent of the 
retrieved images with the largest rms of sky counts, which 
we find is sufficient to reject apparently flawed exposures, 
the images are zero-shifted and scaled according to the sky- 
count statistics and then stacked by the inverse-of-variance 
weighted average using IRAF. The photometric calibration of 
the stacked images is achieved by comparing the stellar pho- 
tometry with that of the reference frames. We find that the 
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Figure 1. Comparisons of stellar magnitudes in the stacked and 
reference frames in u (top left), g (top right), r (middle left), i 
(middle right) and z (bottom left) bands. The dashed lines rep- 
resent the locus where the two measurements are identical. The 
rms errors of the photometric calibration are calculated in the 
magnitude ranges shown by the dotted lines, where the stellar 
photometry is most reliable. 



stellar magnitudes on the stacked and reference frames are 
in excellent agreement, with rms less than 0.05 mag (Fig. 1). 
Fig. 2 shows the comparison of the original and stacked r- 
band images of the same field. More than 100 original SDSS 
frames contribute to each of the stacked frames, and the 
latter images are on average >2mag deeper than the former 
images. 

We run the SOURCE EXTRACTOR on the stacked images to 
extract detected sources. The groups with four or more pix- 
els whose counts are 1.5a above the local background level 
are identified as sources. For every detection, we measure the 
aperture magnitudes in the 4.5-, 3.3-, 2.8- and 2.6-arcsec di- 
ameter apertures as we did for the near-infrared images. The 
photometry errors are calculated from the source photon 
counts and the background noise. The extracted sources are 
cross-identified with the J^-band sources within the maxi- 
mum paring tolerance of 1.0 arcsec. Owing to the deep stack- 
ing of the SDSS images, nearly 90 per cent of the if-band 
sources have counterparts in the g, r, i and z bands. Nearly 




Figure 2. Original (top) and stacked (bottom) SDSS r-band im- 
ages of the same field. 



40 per cent of the if-band sources also have counterparts in 
the u band. 



2.3 Optical Spectroscopy 

We exploit the two redshift surveys carried out on the SDSS 
souther n equatorial stripe; t he VIMOS-VLT Deep Survey 
(VV DS;lLe Fevre et alj[2005h and the DEEP 2 Redshift Sur- 
vey l|Davis et alj|2003tl . Among the four fields of the VVDS 
'Wide' survey, we use the 4-deg 2 field of the F22 (2217+00), 
which lies on the celestial equator, centred at RA 22 h 17 m 
50 s . 4 and Dec. +00° 24 m 00 s . This field coincides with the 
UKIDSS DXS VIMOS 4, and is observed in both the LAS 
and the DXS. We use Field 4 (RA 02 h 30 m , Dec. +00° 00 m ) 
of the DEEP2, one of the T-h survey' fields placed on our 
LAS field. The VVDS adopts a pure /-band flux-limited se- 
lection of the sample while the DEEP2 imposes strict colour 
pre-selection on the spectroscopic targets to favour galax- 
ies at z > 0.7; thus the two surveys are complementary in 
terms of the sample selection. We use only the spectroscopic 
sample with high redshift-quality flags, zflag/zQ — 3 or 4 
for the VVDS/DEEP2. As a result, we obtain 253 LAS K- 
VVDS and 375 LAS A"-DEEP2 galaxies, as well as 1084 
DXS i-f-VVDS galaxies. We note that the redshift surveys 
are deep enough that essentially all the LAS K sources could 
be sampled. 



2.4 Star/Galaxy Classification 

The isT-band sources separate clearly into stars and galaxies 
on the r — z versus z — K diagram as shown in Fig. 3. The 
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Figure 3. The r — z versus z — K diagram for a subset of _R"-band 
sources (black), on which the VVDS stars (blue) and galaxies 
(red) are superimposed. The green line represents the adopted 
star /galaxy classification criteria. 



colours are measured with the 2.8-arcsec diameter aperture 
magnitudes. We define the demarcation between stars and 
galaxies along the minimum surface density on this diagram, 
i.e. the sources redder than z — K = 0.52(r — z) + 1.74 are 
classified as galaxies. The additional criterion of r ~ z < 4 is 
set for galaxies to exclude cool dwarf stars from the sample. 
We obtain 259 082 galaxies with these classification criteria. 
The VVDS classification based on spectra confirms that the 
above scheme works very well, yielding a rate of misclassi- 
fication (stars classified as galaxies and vice versa) of less 
than ~1 per cent. The K-band sources without r- and/or 
z-band detections are excluded from the sample, since their 
extremely red colours suggest that they are mostly galax- 
ies beyond z = 1. Actually, we find that more than 95 per 
cent of the LAS A"- VVDS and LAS A-DEEP2 galaxies are 
detected in both r and z bands in any redshift bins at z ^ 1. 

We show the A-band differential number counts of 
the extracted galaxies in Fig. 4. The y are in excellent 
agreement with p revious measurements (|Daddi et al.ll2000l ; 
iHuang et aT]|200ll ) down to the limiting magnitude of K to t = 
17.9 mag after the detection-completeness correction is ap- 
plied. This suggests that we have successfully constructed a 
well-defined sample of galaxies through the above processes. 



3 REDSHIFT AND STELLAR MASS 
MEASUREMENTS 

3.1 Photometric Redshift 

We estimate the redshifts of galaxies from the observed 
broad-band colours, measured in the 2.8-arcsec diameter 
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Figure 4. The _R"-band differential number counts of the 
extracted galaxies with (squares) and without (diamonds) 
detection-completeness correction. The error bars denote Pois- 
son noise. The blue circles and th e red triangles represent the 
measurements of I Daddi et al.l l|2000l 'l and IHuang et al.l (|200lh , re- 
spectively. 



apertu re, with the optim ized template-fitting method fol- 
lowing |llbert^et^l] l)2006l ). We present a short summary of 
the method below, while the full description of the concept 
and details can be found in the above reference. 

We choose the DXS 7f-VVDS galaxies to optimize the 
spectral templates, leaving the LAS K-VVDS and LAS 
Jf-DEEP2 galaxies as the evaluation sets of the redshift 
measurements (Table 1). First, the DXS A"- VVDS galax- 
ies are classified into five spectral types, Ell, Sbc, Scd, Irr 
and starburst (SB) by least-^ 2 fitting with the appropri- 
ate amounts of dust extinction assuming the exti n ction laws 
of the Small Magellanic Cloud (SMC) by [Pel il 19921) for 
Scd an d Irr and that of starburst galaxies by Calzetti et al.l 
( 2000 ) for SB. The initial spectral templates are taken from 
IColeman. Wu. fc Weedmanl (ll980T) for Ell, Sbc, Scd and Irr 
and from Kinney et al. ( 19961 ) for SB. Then, for each filter 



/, we minimize the sum 



galaxy 

J 



Ax F f , , 

model 



(i) 



where F* bs and cr bs are the observed flux and its error in the 
filter /. The sum is taken over all the sample galaxies. The 
parameters i^ odel and A represent the best-fitting template 
flux and its normalization factor taken from the initial least- 
X 2 fitting. The last term s-^ is a free parameter. While this 
term should be zero in the case of a completely random 
uncertainty in the photometry, we find that it has a non- 
zero value in every filter. These values are at most 0.05 mag 
and are comparable with the expected uncertainty in the 
photometric zero-point calibration. 

The DXS K-VVDS galaxies are re-classified into five 
spectral-type groups with the terms s-^ considered. In each 
group, the observed broad-band fluxes of galaxies are con- 
verted to the rest frame according to the spectroscopic red- 
shifts, after being normalized and de-reddened by the best- 
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Table 1. Summary of the spectroscopic sample. 



Sample 




Number 


Use* 1 


DXS K 


- VVDS 


1084 


training 


LAS K 


- VVDS* 2 


253 


evaluation 


LAS K 


- DEEP2 


375 


evaluation 



* 1 Use in the photometric-redshift measurements. 
* 2 Approximately 70% of the LAS K - VVDS galaxies are also 
the members of the DXS K - VVDS galaxies. 



fitting normalization factor A and the dust extinction. Since 
the galaxies have various redshifts, this conversion gener- 
ates a continuous spectral energy distribution for each spec- 
tral type of galaxies over the relevant rest-frame wavelength 
range. We sort the rest-frame fluxes according to their wave- 
lengths and bin them by groups of points, and connect the 
median flux in each bin to produce the optimized templates. 
We keep the extrapolations provided by the initial templates 
in ultraviolet and infrared wavelengths where no broad-band 
data are available. The starburst template is not optimized, 
in order to retain the emission lines in the template. Finally, 
these optimized templates are interpolated to produce a to- 
tal of 62 templates, the first being Ell and the last being SB, 
to improve the sampling of the redshift-colour space. Below 
we define the spectral type of each galaxy using the best fits 
from among these 62 templates. 

The created spectral templates are fitted to the ob- 
served colours of the actual sample to measure their red- 
shifts. We evaluate the measurement accuracy by applying 
the same procedure to the LAS if -VVDS and LAS if- 
DEEP2 galaxies, as well as the DXS if- VVDS galaxies. The 
DXS if -VVDS galaxies are reduced in number according to 
the LAS if-band detection completeness and are given ad- 
ditional random photometry errors in order to simulate the 
LAS if-band galaxies. We find that the photometric redshifts 
(Zphot) are well correlated with the spectroscopic redshifts 
(zspec) as shown in Fig. 5, thanks to our wide and relatively 
fine wavelength coverage in the u through K bands. The 
deviation between the two (photometric and spectroscopic) 
measurements closely follows a Gaussian distribution with 
standard deviation cr Az /(i+ z ) ~ 0.04 (Az — z Bpcc — z phot ) for 
all three sets of the evaluation sample. Note that ~70 per 
cent of the LAS if -VVDS galaxies are also members of the 
DXS if -VVDS galaxies and account for about a quarter of 
the sample used to build up the spectral templates, so that 
the LAS if -VVDS galaxies do not provide a completely in- 
dependent test of the photometric-redshift accuracy. As a 
further test, we created another template set from the DXS 
if -VVDS galaxies omitting these LAS K -VVDS galaxies 
and repeated the photometric-redshift measurement. This 
test again gives <taz/(i+z) ~ 0.04, which indicates that it 
is a robust estimate of the photometric-redshift uncertainty. 
We show the uncertainty as a function of redshift and stellar 
mass (as determined below) in Table 2. They are relatively 
large in the lowest and highest redshift bins for larger stellar 
mass classes, for which relatively small numbers of sample 
contribute to the spectral templates. We also show the un- 
certainty as a function of spectral type in Table 3, which 
suggests there is little variation of uncertainty among the 
different spectral types. 
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Figure 5. Comparison of the photometric and spectroscopic red- 
shift measurements of LAS if- VVDS (light blue), DXS if- VVDS 
(blue) and LAS i\-DEEP2 (red) galaxies. The solid line shows 
the locus where the two measurements are identical. 



Table 2. Uncertainty (ctaz/(i+z)) m photometric redshift. 









log M* 






Redshift 


10.0 - 10.5 


10.5 - 11.0 


11.0 - 11.5 


11.5 - 12.0 


0.2 - 


0.4 


0.037 (50) 


0.046 (51) 


0.043 ( 6) 


-(0) 


0.4 - 


0.6 


0.040 (35) 


0.037 (55) 


0.039 (25) 


-(o) 


0.6 - 


0.8 


0.021 (35) 


0.030 (73) 


0.032 (93) 


0.034 (20) 


0.8 - 


1.0 


0.033 (17) 


0.038 (71) 


0.046 (55) 


0.047 ( 7) 



Note — Numbers in parentheses represent the number in the 
evaluation samples. 

3.2 Stellar Mass 

The stellar masses (Ai*) of galaxies are determined by fit- 
ting to th e observed colours the stell ar population synthesis 
models of iBruzual fc CharlotJ l|2003h . The aperture magni- 
tudes (m ap ) of the 4.5-, 3.3-, 2.8- and 2.6-arcsec apertures 
are used for the fitting of galaxies with photometric redshifts 
z = 0.2-0.4, 0.4-0.6, 0.6-0.8, and 0.8-1.0, respectively, so 
that we sample the stellar populations consistently within 

Table 3. Uncertainty (o"a z /(i+ 2 )) in photometric redshift. 



Spectral type °"Az/(l+z) 



Ell - 


Sbc 


0.035 (386) 


Sbc - 


Scd 


0.037 (171) 


Scd - 


- Irr 


0.036 ( 35) 


Irr - 


SB 


-(o) 



Note — Numbers in parentheses represent the number in the 
evaluation samples. 
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the central ~20 kpc of all the galaxies. The resultant stel- 
lar mass is then scaled by \Q 0A y K ^p~ Kt ° t ^ to correct for 
aperture loss. We adopt the standard configurations with 
the Padova 1994 stellar evolutionary tracks and B aSel 3.1 
spectral library for the iBruzual fc Charlotl (|2003l ) models. 
We assume three values of metallicity: 0.2, 1 and 2.5 Zq, 
where Zq is the solar metallicity. The star-formation his- 
tory is assumed to take the exponentially declining form 
t -1 exp(— t/r ), where the e-f olding time t is a free param- 
eter, with the ISalpeterl (|l955h initial mass function (IMF). 
Our stellar mass estimates can be approximately converted 
to tho se with another commonly used IMF, that of lChabrierl 
(2003), by adding ~0.25 dex. Other free parameters are the 
age t of the stellar population and the colour excess E(B — V) 
due to the dust extinction of the s tellar radiation, following 
the SMC extinction curve of [Pel (|l992T ). These parameters 
are varied over the plausible ranges of 10 Myr ^ r 10 
Gyr (Alog r = 0.2), 10 Myr sC t sC 10 Gyr (Alog t = 0.1), 
and 0.0 sC E(B - V) < 0.5 mag (AE(B - V) = 0.05), and 
the best-fitting parameter set is searched for by the least-x 2 
method for each galaxy (the values in parentheses represent 
the grid intervals). The additional error of 0.05 mag is added 
in quadrature to all band magnitudes in the fitting in order 
to take into account the uncertainty in the photometric zero- 
point calibration. 



We derive two kinds of stellar mass for the spectro- 
scopic sample, i.e. the stellar mass with spectroscopic red- 
shifts (Ad*, spec) and the stellar mass with photometric red- 
shifts (M+^hot). The difference between the two measures, 
AlogM* = logM* !S pe C — logA/*, p hoi, is found to be clearly 
correlated with the photometric redshift deviation Az = 
Zspcc — Zphot- Such a correlation is expected, since a larger 
Zphot leads to a larger estimate of the galaxy luminosity, 
which then leads to a larger estimate of stellar mass M* iP h ot . 
The observed relation between Az and AlogM* is actually 
quite consistent with this expected correlation. Another ex- 
pected cause of the correlation between Az and AlogA/* 
comes from the fact that larger estimates of z p hot lead to 
systematically shorter rest-frame wavelengths to which each 
of the observing wavebands corresponds. After removal of 
the above first component of the systematic correlation, we 
found marginal evidence for the second correlation in our 
sample, which is AlogAf* = +/- (0.04 ± 0.07) when Az is 
negative/positive. In addition, we consider the uncertainty 
associated with the least-x 2 model fitting. It is evaluated 
by the la confidence surface of the \ 2 distributions in the 
fitting parameter space. 



We show the total amplitudes of the stellar mass uncer- 
tainty ((TAiogA/*) as a function of redshift and stellar mass 
in Table 4. Those as a function of spectral type are shown 
in Table 5. The above estimates of error amplitudes and the 
correlations between Az and AlogM* are taken into account 
in the Monte Carlo simulation presented below. Note that 
further different assumptions on the stellar mass estimation, 
such as different stellar population synthesis models and dif- 
ferent IMF, can cause additional uncertainty in the derived 
properties of galaxies. We will address this issue in Section 
5. 



Table 4. Uncertainty (ffAlogM, ) i n stellar mass. 









log Ah 






Redshift 


10.0 - 10.5 


10.5 - 11.0 


11.0 - 11.5 


11.5 - 12.0 


0.2 - 


0.4 


0.19 (50) 


0.20 (51) 


0.22 ( 6) 


-(0) 


0.4 - 


0.6 


0.19 (35) 


0.17 (55) 


0.16 (25) 


-(0) 


0.6 - 


0.8 


0.12 (35) 


0.16 (73) 


0.17 (93) 


0.19 (20) 


0.8 - 


1.0 


0.14 (17) 


0.17 (71) 


0.19 (55) 


0.25 ( 7) 



Note — Numbers in parentheses represent the number in the 
evaluation samples. 

Table 5. Uncertainty (oAiogM* ) i n stellar mass 



Spectral type o"AiogM* 



Ell - 


Sbc 


0.16 (386) 


Sbc - 


Scd 


0.17 (171) 


Scd - 


- Irr 


0.17 ( 35) 


Irr - 


SB 


-(o) 



Note — Numbers in parentheses represent the number in the 
evaluation samples. 

4 RESULTS 



We define our massive galaxy sample using two stellar mass 
classes, i.e. 1O 11O ~ 11J5 M galaxies with 10 110 M Q < M* < 
10 115 A/ Q and 10 11 - 5 " 120 M Q galaxies with IO^^Mq < 
AL < 1O 12O M . The galaxies are grouped into four red- 
shift bins, z = 0.2 - 0.4, 0.4 - 0.6, 0.6 - 0.8, and 0.8 - 1.0. 
Total numbers included in the sample are summarized in 
Table 6. The median photometry errors in the r-, z-, and 
A'-band aperture magnitudes and in the if-band total mag- 



nitudes (r a 



Zap , -^"ap , -^tot ) 



are also listed. The aperture 



magnitude errors are generally smaller than the typical un- 
certainty in the photometric zero-point calibration (~0.05 
mag). 

We show the differential number counts of the massive 
galaxies in Fig. 6. It shows that the number-count distribu- 
tions of the 10 11,5-12 '°Mq galaxies have faint-end drop-offs 
at magnitudes brighter than the limiting magnitude, which 
assures the near-complete detection of this population. On 
the other hand, the faintest of the W 11 '°~ 115 Mq galaxies at 
high redshifts (z > 0.6) fall below the limiting magnitude, 
and are thus left uncounted. In order to estimate the lost 
fraction of 1O 11O_11 ' 5 M0 galaxies at these redshifts, we de- 
rive the detection completeness specifically for these galax- 
ies as follows. We take each of the 1O 11O_11 ' 5 M0 galaxies 
in the lowest redshift bin and assign random redshifts in 
the 0.6 < z < 0.8 and 0.8 < z < 1.0 ranges. The galaxies 
are dimmed and reduced in apparent size according to the 
assigned redshifts, placed on random positions of the LAS 
K-band images, and then extracted by SOURCE EXTRACTOR in 
the same way as the actual sample sources are detected. The 
recovery rate of the embedded objects as a function of their 
magnitudes gives the detection completeness of the galaxies, 
which we find is significantly better than that of the whole 
sample derived before. The 50 per cent detection complete- 
ness is actually achieved at K to t = 18.6 mag instead of K to t 
= 17.9 mag, and almost all galaxies brighter than K to t = 
17.9 mag are detected. With the new detection-completeness 
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Table 6. Summary of the massive galaxy sample. 







log Mi, = 


11.0 - 11.5 


log Mi, = 


11.5 - 12.0 


Redshift 


Number 


Photometry Error* 


Number 


Photometry Error* 


0.2 - 


0.4 


9,720 


(<0.01, <0.01, 0.01, 0.05) 


1,408 


(<0.01, <0.01, <0.01, 0.03) 


0.4 - 


0.6 


15,300 


(<0.01, <0.01, 0.02, 0.08) 


572 


(<0.01, <0.01, 0.01, 0.04) 


0.6 - 


0.8 


18,582 


(0.02, 0.02, 0.03, 0.13) 


815 


(0.01, 0.01, 0.02, 0.09) 


0.8 - 


1.0 


12,371 


(0.04, 0.02, 0.04, 0.16) 


613 


(0.03, 0.02, 0.03, 0.12) 



Note (*) — Median photometry errors in the r-, z- and K-band aperture magnitudes and in the K-band total magnitudes (r ap , Zap, 

K&p, Ktot)- 



< 0.6 




14 15 16 17 1 
K (mag) 



14 15 16 17 18 19 
K (mag) 



Figure 6. The differential number counts of 10 110 115 Mq 
(green) and 1O 11 ' 5_12 O M0 (orange) galaxies at z = 0.2-0.4 (top 
left), 0.4-0.6 (top right), 0.6-0.8 (bottom left) and 0.8-1.0 (bot- 
tom right). The black dashed lines represent our formal limiting 
magnitude of K tot = 17.9 mag. 



function taken into account, the fractions of 1O 11O_11 ' 5 M0 
galaxies fainter than the formal limiting magnitude (K = 
17.9 mag, thus uncounted) are 11 and 16 per cent at z = 
0.6-0.8 and 0.8-1.0, respectively 

We estimate the uncertainty in the measured numbers 
of massive galaxies by a Monte Carlo simulation, as follows. 
First, we generate a mock galaxy catalogue containing 10 
galaxies for each stellar mass (AlogA/* = 0.1) and redshift 
(Az = 0.02) bin in the ranges 8.0 < log M* < 13.0 and 
0.0 < z < 1.4, where the numbers in parentheses represent 
the bin widths. Each galaxy is assigned a weighting factor 
corresponding to the num ber density of gal axies, following 
the galaxy mass function of lCole et al.l |2001). Next the red- 
shifts of galaxies are given perturbations following the mea- 
sured uncertainty of the actual sample. The stellar masses 
are also given perturbations, correlated with the photomet- 
ric redshift perturbations as explored in Section 3.2. Then 
the mock galaxies are weighted by their weighting factors 
and redistributed, and counted in the redshift and stellar 



mass bins to obtain the output mass function. We repeat 
the calculation 100 times, varying the random components. 

Fig. 7 shows the results of the simulation. We observe 
both systematic and random components in the resultant 
error estimates. The systematic component is evident in 
Mj, > 10 12 ■°M fl bins. This is the so-called Eddington bias 
(Eddi ngtonl Il913l ). caused by the steep slope of the high 
end of the mass function; simply put, a small portion of 
the less massive, much more numerous galaxies could con- 
taminate the more massive classes owing to measurement 
errors, which significantly alters the steep part of the mass 
function. This is why we limit our massive galaxy sample 
to those with M 4 < 1O 12 '°M0; the systematic increases in 
number are found to be insignificant for our mass ranges, i.e. 
negligible for the W 110 ~ 11 ' 5 Mq galaxies and 40-60 per cent 
for the 1O 11 ' 5_12 O M0 galaxies. The output number densi- 
ties from the 100 repeated calculations scatter around the 
systematic components, which yields the random compo- 
nents of the measurement error. The standard deviations of 
the scatter are ~5 per cent and ~12 per cent of the num- 
bers of 1O 11O_11 ' 5 M and 1O 11 - 5 ~ 12 ' O M galaxies, respec- 
tively. Note that the current estimate is not perfect, since 
we assume a non-evolving galaxy stellar mass function at 
< z < 1, while our results show clear signs of its evolu- 
tion (see below). We will discuss this issue further in the 
following section. 

Another source of uncertainty comes from cosmic vari- 
ance. We estimate cosmic variance by dividing our sample 
into five subfields along the right ascension, ~11 deg 2 each, 
and then calculating the fractional variation of the measured 
numbers of massive galaxies in these subfields. We find that 
the fractional variations are <10 per cent and <20 per cent 
for the 10 n - 0-u - 5 Af s and 1O 115 " 12 O M galaxies, respec- 
tively, in each of the four redshift bins. Considering that the 
total fields are five times larger than the subfields, we con- 
clude that cosmic variance could affect the measured num- 
ber of massive galaxies in each redshift bin by up to 5 per 
cent and 10 per cent for the two classes of galaxies, respec- 
tively. The above estimates are ro ughly consistent with th e 
theoretical predictions provided bv lSomerville et al.l (|2004h . 

We show our results with regard to the number-density 
measurements in Table 7, and plo t them along wit h the mea- 
surements for the local Universe l|Cole et al.ll200ll ) in Fig. 8. 
The error bars take into account all uncertainties considered 
above, as well as the Poisson noises. The number densities 
of the 10 11 - S_12 -°M Q galaxies are corrected for Eddington 
bias, although this correction has little significance for our 
final conclusions. The local densi ties were no r maliz ed to take 
into account the assumptions of ICole et alj (|200ll ) with re- 
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0.2 < z < 0.4 0.4 < z < 0.6 




8 9 10 11 12 13 8 9 10 11 12 13 

log M log M 



Figure 7. The input (black line) and output (100 red lines) 
galaxy mass functions of the Monte Carlo simulation at z = 0.2- 
0.4 (top left), 0.4-0.6 (top right), 0.6-0.8 (bottom left) and 0.8-1.0 
(bottom right). The dashed lines show the mass ranges in which 
our massive galaxy sample is defined. 



Table 7. Number densities of the massive galaxies. 







log M* 




Redshift 


11.0 - 11.5 


11.5 - 12.0 


0.2 


- 0.4 


10.9 ± 0.8 


0.98 ± 0.15 


0.4 


- 0.6 


7.8 ± 0.6 


0.19 ± 0.03 


0.6 


- 0.8 


6.8 ± 0.5 


0.17 ± 0.03 


0.8 


- 1.0 


3.7 ± 0.3 


0.11 ± 0.02 



Note — Number densities are given in units of 10 4 Mpc 3 log 

gard to stellar population synthesis that differ from ours; 
the major difference is that they assume a constant forma- 
tion redshift of galaxies at Zf = 20, while we vary the age of 
the stellar population as a free parameter. The normaliza- 
tion factor, ~ 0.4, is derived by applying their assumption 
to our 1O 11O ~ 11 ' 5 M0 galaxies in the lowest redshift bin z = 
0.2-0.4. As seen in the figure, we find that the most massive 
(1O 11 ' 5_12 O M0) galaxies have experienced rapid evolution 
in number since z — 1. On the other hand, the number den- 
sities of the less massive (10 11 ' 0- M ) systems show a 
rather mild evolution during the same period. 



5 DISCUSSION 

The measured number density evolution of massive galaxies 
shows clear signs of the hierarchical evolution of these sys- 
tems. Such a galaxy evolution scenario is predicted in the lat- 
est galaxy formation models based on ACDM theory. In Fig. 



Redshift 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 
i 1 1 1 1 1 1 — i — i — i — n - ! 




J I I I I I I I I I I I I I I I L 



14 12 10 8 6 

Age of the universe (Gyr) 

Figure 8. The number densities of massive galaxies versus the 
age of the Universe or redshift. The green and orange symbols and 
lines represent the 10 110-11 - 5 M and 10 11 - 5_12 M Q galaxies, 
respectively. The filled circles represent our measurements while 
the open circles represent measurements in the local Universe 
l|Cole et al.ll20oih . The solid lines represent the predictions of the 
Millennium Simulation with the semi-analytic galaxy formation 
model (see text). 

8 we overlay the predictions of the Millennium Simulation 
(Lemson et al. 2006), the largest numerical simulation to 
date based on AC DM theory, with the s emian alytic galaxy 
formation model of iDe Lucia fc Blaizotl (|2007h . scaled to fit 
to the local observations (by a factor of 0.4) . The stellar mass 
of the Millennium model has been shifted by +0.25 dex in 
order to c orrect for the di fferent IMFs adopted (the model 
adopts the lChabrierl l|2003l ) IMF). The observed hierarchical 
pattern of evolution is consistent with the prediction of the 
model, while we find some discrepancies between the obser- 
vation and the model (e.g. the 1O 11O_11 ' 5 M0 galaxies at z — 
0.8-1.0). This indicates that the basic idea of the bottom-up 
construction of galaxy systems is valid at least for the most 
massive galaxies with M* > 1O 11 M0. 

Meanwhile, we point out that our measurements plot- 
ted in Fig. 8 seem to have an apparently unnatural feature 
at z = 0.2-0.4, where the number densities are significantly 
high relative to the overall trend considering their estimated 
errors, implying residual systematic effects from the Edding- 
ton bias. We check this by investigating the systematic in- 
creases in the number of massive galaxies in the spectro- 
scopic sample, from those obtained with the spectroscopic 
redshifts to those with photometric redshifts (we consider 
only those redshift and stellar mass classes with more than 
five spectroscopic sources). The above Monte Carlo simu- 
lation shows that the use of photometric redshifts (and the 
consequent stellar mass fluctuations) is the dominant source 
of the Eddington bias. We list the measured systematic in- 
creases in Table 8 separately for the VVDS and DEEP2 
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galaxies, since the amplitudes of the Eddington bias are 
subject to the redshift distribution of sources. The DXS 
7f-VVDS galaxies have been given additional photometry 
errors and detection incompleteness in order to simulate the 
LAS K galaxies. The systematic increases estimated in the 
Monte Carlo simulation are also listed in Table 8. The table 
shows that the systematic increases measured in the spec- 
troscopic sample are mostly less than the estimates from 
the Monte Carlo simulation, while the one for the VVDS 
10 11-0_11 - 6 Mq galaxies at z = 0.2-0.4 is exceptionally high. 
While it is based on a small sample, it provides marginal ev- 
idence that the Eddington bias is unexpectedly significant 
in this lowest redshift bin, due to some unquantified sources 
of systematic uncertainty. In that case, the number densities 
measured at z = 0.2-0.4 should be regarded as the upper 
limits. 

We also note that we adopt a non-evolving galaxy stel- 
lar mass function in the Monte Carlo simulation, while our 
results suggest the steepening of the high end of the mass 
function toward high redshifts; thus the estimated amount 
of Eddington bias should be regarded as a lower limit. When 
such a steepening of the mass function is taken into account 
in the correction of the Eddington bias, we obtain even lower 
numbers of galaxies in the more massive classes at higher 
redshifts, whi ch would further s t rengt hen our conclusion. 

Recently, iMarchesini et all (|2009l ) provided a detailed 
analysis of random and systematic uncertainties affecting 
the galaxy stellar mass function. They adopt 14 differ- 
ent stellar mass estimations with different combinations of 
metallicity, dust extinction law, stellar population synthesis 
models and IMF, and show that the derived number den- 
sity of galaxies in a given stellar mass bin could be altered 
by up to 1 dex. The 'bottom-light' IMFs in partic ular, with 
a defi cit of low-mass stars relative to a standard IChabrieil 
|2003t ) IMF, give significantly different results for the stellar 
mass function from other classical IMFs. A more top-heavy 
IMF at hi g her re dshif ts is actually s ugges ted by, for exam- 
ple, iDavel 4 2008) and Ivan Dokkuml |200ct ). Thus we point 
out that our results are subject to a systematic change of 
these stellar-population properties during < z < 1. Future 
improvements in stellar population models based on new ob- 
servations are eagerly awaited to overcome these large un- 
certainties inherent in stellar mass measurements. 

In order to probe the star-forming properties of massive 
galaxies, here we investigate their rest-frame optical colours. 
Nearby galaxies are known to show a clear bimodality in 
the optical colour distribution, in which early-type galax- 
ies form a narrow red sequence that is separated from blue 
star-forming populations by a valley of the galax y distri- 
bution (e.g-. IStrateva et al.lHoOll : Irlogg et alj|2003h. A sim- 
ilar bimodality is observed out to z > 1 (e.g.. iLin et al.l 
Il999l : llm et"al]|2002l ; lBell et al.ll2004h . We calculate the rest- 
frame U- and U-band magnitudes of massive galaxies by k- 
correcting the nearest observed r, i or z-band magnitudes, 
where the amounts of the fc-corrections are estimated from 
the best-fitting spectral population synthesis models. We 
show the resultant colour distributions in Fig. 9. 

As Fig. 9 shows, the less massive (< 1O 11O M ) galaxies 
show a clear colour bimodality, as expected, with the peak 
colour of the red sequence at U — V ~ 1.2 and a valley 
of galaxy distributions at U — V ~ 1.0 in all redshift bins. 
Compared with these galaxies, massive (> 10 °M ) galax- 
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Figure 9. U — V colour distributions of the massive galaxies at z 
= 0.2-0.4 (top left), 0.4-0.6 (top right), 0.6-0.8 (bottom left) and 
0.8—1.0 (bottom right). The grey, green and orange lines represent 
the < lO n - o Af , 10 110 - n - 5 M Q and lO n - 5 - 12O M galaxies, 
respectively. The dotted lines show the demarcation between the 
blue and red populations. 



ies are apparently dominated by the red population, with 
conspicuous peaks at U — V ~ 1.2. We divide the galaxies 
into blue and red populations at U — V — 1.0 and calculate 
the blue fractions (fractions of the blue population). The 
associated errors are estimated by repeatedly giving ran- 
dom fluctuations to the U — V colours, taking into account 
the uncertainties in the photometry and the amount of k- 
correction, and then re-measuring the blue fractions. We find 
that the fluctuated blue fractions are systematically larger 
than the original values, since there is a greater red popu- 
lation than blue around the demarcation U — V = 1.0, and 
correct for the effect (the correction amounts are included 
in the final errors). We plot the measured blue fractions as 
a function of redshift and stellar mass in Fig. 10, and also 
list them in Table 9. One can see that the blue fractions are 
significantly lower in more massive galaxies, and that the 
fractions in massive systems decrease toward the local Uni- 
verse. The blue fractions in < lO llo Af galaxies increase 
from z = 0.4-0.6 to 0.2-0.4 because the dominant popula- 
tion within the sample shifts to bluer, less massive galaxies 
toward the local Universe due to the fainter detection limit. 
In fact we observe a decreasing blue fraction toward the local 
Universe, as seen in massive galaxies, if we take the subsam- 
ple with stellar mass 1O 1O - 5 M < M* < 1O 11O M (see Table 
9). As discussed above, the 1O 11O_11 ' 5 M galaxies at z — 
0.2-0.4 and the 10 115 " 12 A'/ galaxies in all redshift bins 
could have considerable fractions of contamination from less 
massive galaxies, which likely have bluer U — V colours. Ac- 
tually, investigating the spectroscopic sample shows that the 
contamination makes the mean U — V colours bluer by up 
to 0.1 mag. Thus the true blue fractions in the above classes 
of galaxies could be even smaller than the present measure- 
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Table 8. Estimates of the Eddington bias* 1 . 



log Af* = 11.0 - 11.5 log Mi, = 11.5 - 12.0 



Redshift 


MCS* 2 


VVDS* 3 


DEEP2* 3 


MCS* 2 


VVDS* 3 


DEEP2* 3 


0.2 - 0.4 


< 0.1 


0.7 (6) 


-(0) 


0.6 


- (0) 


- (0) 


0.4 - 0.6 


< 0.1 


< 0.1 (23) 


- (2) 


0.5 


- (0) 


- (0) 


0.6 - 0.8 


0.1 


0.1 (21) 


< 0.1 (72) 


0.6 


0.3 (7) 


0.3 (13) 


0.8 - 1.0 


< 0.1 


- (2) 


< 0.1 (53) 


0.4 


- (0) 


0.1 (7) 



^Rates of increase (increased amounts divided by the original values) of the number densities are listed. 
* 2 Estimates from the Monte Carlo simulation. 

* 3 Estimates from the VVDS and DEEP2 spectroscopic samples (the number in the samples is shown in the parentheses). 



Table 9. The blue population fraction. 









log M+ 




Redshift 


< 11.0 (10.5 - 11.0) 


11.0 - 11.5 


11.5 - 12.0 


0.2 - 


0.4 


0.54 ± 0.01 (0.36 ± 0.01) 


0.12 ± 0.05 


< 0.06 


0.4 - 


0.6 


0.47 ± 0.01 (0.43 ± 0.01) 


0.25 ± 0.03 


0.18 ± 0.03 


0.6 - 


0.8 


0.51 ± 0.01 (0.43 ± 0.01) 


0.28 ± 0.02 


0.19 ± 0.06 


0.8 - 


1.0 


0.59 ± 0.01 (0.56 ± 0.01) 


0.29 ± 0.02 


0.21 ± 0.03 



merits. The lower blue fractions in more massive galaxies and 
the decreasing trend toward the local Universe implies ma- 
jor star formation at higher redshifts, which is in line wit h 
'downsizing' of the star formation (e.g.. lCowie et al1ll996f ). 

The above measurements suggest that the majority 
of massive galaxies are fairly quiescent, while of the rest 
a considerable fraction of galaxies are experiencing ac- 
tive star formation, especially at higher redshifts (~30 per 
cent at z ~ 1). Such activ e star formatio n in m assive 
galaxies is also reported in IConselice et al.l (|2007| ). who 
find that nearly half of their massive (M* > 1O 11 M0) 
galaxies at 0.4 < z < 1.4 are detected in the Spitzer 
Space Telescope/MlPS 24-/im band and the average star- 
formation rate amounts to ~ 50Mq yr _1 . Our measure- 
ments of blue fractions indicate that the star-formation 
activity in massive galaxies is gradually quenched toward 
the local Universe, leaving the most massive galaxies on 
the red sequence. Star-formation quenching processes above 
a certain stellar mass limit are actually proposed, such 
as the internal feedback of mass assembly caused by ac- 
tive g alactic nuclei (e.g., Silk fc Reesl 19981; Granato et al.l 



12004 ISpringel. Pi Matteo. fc HernquistJ I2005T I. The above 
scenario is consistent with our primary results that the ac- 
tive bottom-up formation of massive galaxies is going on 
during < z < 1. 

Finally, we comment on the compatibility of the present 
results with previous studies. There are a num ber of studies 
of luminou s red galaxies (LR Gs) at z < 1 (e.g.. lBrown et all 
l200l 120081 : ICool et al.ll2008h covering up to ~10 deg 2 . Au- 
thors consistently suggest that the LRGs show little evolu- 
tion in number density since z ~ 1. However, the mass-to- 
optical luminosity ratio of galaxies has a significant scatter 
even for the massive systems, so that galaxies with a cer- 
tain stellar mass are not quite the equivalent population 
of galaxies with a certain optical luminosity. This leads to a 
consequence of most significance for the most massive galax- 
ies: a small portion of the less massive, much more numerous 
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Figure 10. The fraction of the blue population in < 1O 11O M0 
(grey solid), 10 11 - 0-11 ' 5 A/ Q (green) and 10 11 - 5-120 M (orange) 
galaxies, respectively, as a function of redshift. 



galaxies could contaminate the luminous class of galaxies if 
their mass-to-luminosity ratios were slightly less than the 
average, and thus could easily dominate the luminous pop- 
ulation. Therefore a subtle (in absolute amplitude) change 
in the number of most massive galaxies could be drowned 
out in the measured evolution of the LRGs. 

Studies also exist of massive galaxies at z < 1 with 
stellar mass measurements based on infrared photome- 
try (e.g.. lConselice et alj|2007l ; lllbert et al.ll2009l ). although 
these studies cover a much smaller field of view (< 1.5 deg 2 ) 
than ours. In contrast to the present results, they report lit- 
tle evolution in number of the most massive (> 10 11 ' 5 Mq) 
galaxies. At least part of the discrepancy could be due to 
small-number statistics and cosmic variance. Actually, while 
we find ~1000 samples of the most massive galaxies in each 
redshift bin from our 55.2 deg 2 , the number of samples ob- 
served over ~1.5 deg 2 should be only ~30. We estimate the 
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effects of cosmic variance by dividing our total field into 
small subfields, each covering ~ 1.5 deg 2 , and measure the 
number-density fluctuations of massive galaxies among the 
subfields. As a result, we find that the number densities of 
the most massive galaxies measured over ~ 1.5 deg 2 can fluc- 
tuate by up to a factor of a few. However, we are not sure 
whether the above uncertainties alone can account for the 
discrepancy between the present results and previous ones. 
The number-density measurements at the steep high end of 
the galaxy stellar mass function could be heavily affected by 
contamination arising from less massive galaxies, thus quite 
accurate analysis is required to unveil the subtle evolution 
of the most massive galaxies. 

In essence, our measurements provide a unique oppor- 
tunity to investigate the mean properties and evolution of 
the most massive galaxies, owing to the reliable estimates 
of photometric redshifts and stellar masses conducted over 
an unprecedentedly large field of view. What is observation- 
ally clear is that we have discovered a substantial deficit of 
the most massive galaxies out to z = 1 compared with the 
local Universe. The analysis of the rest-frame U — V colour 
distributions indicates that star-formation activity might be 
responsible for the active build-up of these systems, while it 
is possible that a s o-called dry merger (e.g.. lBell et aljfeooj ; 
van Dokkuml boOS) is the main driver of the evolution. Ac- 
tually, some observations suggest that a substantial fraction 
of massive early-type galaxies go through active evolution 
in terms of the galaxy structure as well as the star forma- 
tion since z ~ 1 (e.g.. iTreu et al.ll2005l ; Ivan der Wei et al.l 
2008). The present results provide crucial evidence of hier- 
archical galaxy formation, the missing piece of observation 
required to chart a course for future theoretical models based 
on ACDM theory. 



6 SUMMARY 

We present an analysis of ~60 000 massive galaxies with stel- 
lar masses 10 n M Q < K'L < 10 12 M Q in an unprecedentedly 
large field of view of 55.2 deg 2 . The galaxies are drawn from 
the UKIDSS Large Area Survey K-band images on the SDSS 
southern equatorial stripe. We have created deep-stacked u, 
g, r, i and z-band images from the SDSS Supplemental and 
Supernova Survey image frames, which results in ~90 per 
cent counterparts of the Tf-band sources in the g, r, i and 
z bands. We also exploit the redshift surveys conducted on 
the SDSS southern equatorial stripe, namely the VIMOS- 
VLT Deep Survey and the DEEP2 Redshift Survey, in or- 
der to obtain accurate photometric redshifts and associated 
uncertainties for the galaxies. Stellar masses are estimated 
by comparing the observed broad-band colours with stellar 
population synthesis models. 

In each of the redshift bins z = 0.2-0.4, 0.4-0.6, 0.6- 
0.8 and 0.8-1.0, we obtain ~10 000 and ~1 000 galax- 
ies with stellar masses 10 110 Af e < KL < lO n - 5 M and 
10 115 M Q < M* < 10 12 M Q , respectively. The galaxies are 
almost completely detected out to z = 1, and form by far 
the largest sample of massive galaxies reaching to the Uni- 
verse at about half its present age. We find that the most 
massive (1O 115 M < K'h < 1O 12 O M ) galaxies have ex- 
perienced rapid growth in number since z = 1, while the 
number densities of less massive systems show rather mild 



evolution. Such a hierarchical trend of evolution is consis- 
tent with the predictions of the current semi-analytic galaxy 
formation model based on ACDM theory. While the major- 
ity of the massive galaxies are red-sequence populations, we 
find that a considerable fraction are blue star-forming galax- 
ies. The blue fraction is less in more massive systems and 
decreases toward the local Universe, leaving the red, most 
massive galaxies at low redshifts, which further supports the 
idea of active bottom-up formation of these populations dur- 
ing < z < 1. The present results provide strong evidence 
that galaxy formation proceeds in a hierarchical way, and 
place stringent observational constraints on future theoreti- 
cal models. 
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